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Abstract. The paper presents numerical simulations performed on dielectric 
properties of two-dimensional binary composites on eleven regular space filling 
tessellations. First, significant contributions of different parameters, which play 
an important role in the electrical properties of the composite, are introduced 
both for designing and analyzing material mixtures. Later, influence of structural 
differences and intrinsic electrical properties of constituents on the composite's over 
all electrical properties are investigated. The structural differences are resolved by 
the spectral density representation approach. At low concentrations of inclusions 
(concentrations lower than the percolation threshold), the spectral density functions 
are delta-sequences, which corresponds to the predictions of the general Maxwell- 
Garnctt mixture formula. At high concentrations of inclusions (close to the percolation 
threshold) systems exhibit non-Debye type dielectric dispersions, and the spectral 
density function differ from each other and that predicted by the Maxwell-Garnett 
expression. The analysis of the dielectric dispersions with an empirical formula 
also illustrate structural differences between the considered geometries, however, 
the information is not qualitative. The empirical formula can only be used to 
compare structures. The spectral representation on the other hand is a concrete 
method to characterize the structures of dielectric mixtures. Therefore, like in other 
spectroscopic techniques, a look-up table might be usefull to classify/characterize 
composite materials structure. This can be achieved by generating dielectric data for 
known structures as presented and emphasized in this study. Finally, the numerical 
technique, without any a-priori assumptions, for extracting the spectral density 
function is also presented. 



1. Introduction 

Composite materials are extensively used for various mechanical, thermal and electrical 
application nowadays. Their electrical properties gain significant interest due to their 
possible applications as special materials, sensors and specially due to understanding the 
nature of already existing materials by utilizing electrical impedance measurements 0, 
0, 0|, such as in polymer physics and geophysics. Example of composites can be 
classified as man-made and natural mixtures |4j. In the former class are composites 
(fiber, particulated and percolating), porous materials (gels, cellular solids and foams), 
colloids, microemulsions, block copolymers, concrete and even biological samples, which 
are some sort of emulsion in a continuous background of solution. Some of the natural 



mixtures are polycrystals, soil, sandstone, granular media, and biological materials such 
as wood, bone, blood, lungs, tissues, tumors etc. 

Although bulk structure of materials are not tailored yet, by the introduction of 



micro-electronics and nanotechnology 



0,3, 



we are now able to manufacture desired 



structures using the lithography techniques. It is not yet possible to produce materials 
with three-dimensional structures, however, besides the micro-electronic mechanical 
systems (known as MEMS) attempts have been made to produce materials as sensors 
with lasers and photo-curable solutions 0, 0. The process to miniaturize man-made 
equipment pushes scientist and engineers to come up with new manufacturing ideas. 
Currently most of the research has focused on the surface manipulation of materials, 
however, the electro-kinetic SMEHElEl and acoustic techniques would infact create 
the opportunity to tailor the bulk properties and structures of material composites. 
If^this^ is possible, then computer simulations |bJ, jl4} . 1£ , Eil EU EB, 21, 22, 
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28, 



formulas 





and dielectric mixture 



would be used to predict the behavior 



of materials before and to characterize after the manufacturing. 

The analytical or empirical formulas are general expressions for the mixtures, they 
are only useful, would yield valuable information, when the composite lacks complexity 
in its topology and one of the phases has low volume fractions |24l l25l l26l I32L 145 1 . In the 
last couple of decades with the improvements in the numerical codes and computation 
spe ed, numerical solutions to obtain the electrical properties have been widely used La 
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The numerical calculations overcome the disadvantages of the previously mentioned 
items, however, they have also handicaps such as the size of the problem, the ratio of 
the largest to smallest unit in the structure, and numerical errors. 

One of the novel features in composites, unfortunately not used widely, was 
proposed by Bergman 0, 0, Esl ]. which is called the spectral representation theory. 
In this theory if the electrical properties of the constituents of a binary mixture 
are known, then the contribution of the topology, or in other words the geometrical 
arrangement of constituents, can be extracted from the information hidden in the 
Maxwell- Wagner-Sillars polarization 0, 0, EJ. In addition one can used this 
approach to resolve the dielectric permittivity of a unknown material in a host medium, 
if a test has been performed with a known material with the same geometrical shape 
and distribution as the unknown one in the host medium. This can be achieved by 
using the same spectral density function for the unknown material and reverse engineer 
to obtain the dielectric properties of the unknown material. Previously, researcher 
have proposed cumbersome analytical expressions, which are based on the adopted 



spectral density functions 
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Here, it is not only discussed how to extract the 



unknown spectral density function but it is applied to extract structural differences 
of two-dimensional structures, which are unidirectional fibers in a host medium in 
three-dimensions. The extraction is based on a numerical technique developed by the 



present author 54, 55]. It is based on the Monte Carlo integration hypothesis and the 



constrained-least-squares algorithm. There are no a-priori assumptions in the procedure 
as the previous researchers proposed. The method is tested on eleven two-dimensional 
space-filling tessellation lattices/structures [5 fi^ The dielectric permittivity of structures 
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calculated by the finite element method 0, 
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of the intrinsic electrical properties of the constituents on the dielectric relaxations are 
also investigated. 



2. Background 

2.1. Factors influencing composite properties 

Composites are new sort of materials in which we have the possibility to tailor 
their properties for desired applications. The tailoring can be as in the form of (i) 
combining the intrinsic properties of the constituents, (ii) creating a new non-existing 
property, such as, interfacial polarization in mixtures, electro-mechanical activity in 
porous materials, meta-materials (negative permittivity/permeability and Poisson's 
ratio materials), and (iii) nullifying a property, this is achieved by destroy the structure 
of the continuous matrix. 

If we focus on the significant items influencing the dielectric properties of 
composites, they can be briefly listed as follows; 

• intrinsic properties of constituents, i.e. ohmic conductivity a, high frequency 
dielectric permittivity e and dielectric susceptibility x( u )- The first two quantities 
are material constants, the latter quantity contains the dynamic properties of the 
constituents. The complex dielectric permittivity of a material e can in general be 
expressed as a function of the angular frequency uj 

e{u) = e + x{oj) + aiiSoijy 1 (1) 

where the over-lines represent the quantities are complex numbers. The permittivity 
is a function of angular frequency u>, and losses due to ohmic condiutivity a 
contribute to the imaginary part of the permittivity. Eq is the permittivity of the 
free space, 6q = 8.854 Fm" 1 and i = \J — 1. Keep in mind that the permittivity of 
materials are quantized meaning that the permittivity axis is not continuous. This 
is important while designing high and low permittivity materials due to selection 
of constituents. The conductivity of materials on the other hand can be altered by 
doping. 

• concentrations or volume fractions of constituents. This case is widely used in 
the industry to make semi-conducting/conducting polymers by using carbon-black 
powder, which is a dopant. 

• size, shape, orientation and distribution of inclusions (if one of the phases is in 
particulated form). Influence of particle size on electrical properties have been 
investigate in percolating systems. Although it is advantageous to use mono- 
dispersed spherical inclusions to design composites, the mechanical compatibility is 



reduced due to adhesion. Therefore, usually fillers with arbitrary shapes and largest 
surface area are selected. Still, there is no concrete examples of bulk materials, in 
which orientation and distribution of the inclusions are controlled. 

• mixture topology. There is no concrete examples of controling the geometrical 
arrangements of constituents other than micro-electronics designs for lateral 
structures and arrangement of nano-scale particles on surfaces. The current study 
is focused on this aspect of mixtures, by performing numerical simulations on "ideal 
composites", in which their electrical and structural properties are manupulated. 
When the inclusion phase is expensive control of the distribution of them can in 
principle reduce materials cost or toxic vast without affecting the over-all composite 
properties. 

• inter-phase/interface between the phases. The interface between phases is an 
interesting topic since it is not clear how the transition from one phase to the other 
takes place. One of the open questions is the presence of an inter-phase and its 
dimensions and its influence on the overall composite properties. For example, 
samples prepared with surface modified fillers exhibit changes in the dielectric 
dispersions |5j[|. Moreover, as the particles or phases agglomerate and the sizes 
of the phases decrease the surface to volume ratio of inclusion phase increase and 
the system is dominated with the surface effect. 



2.2. Importance of the mixture topology 

As an example, if we are asked to come up with a novel materials design for a specific 
application, we must ponder about the items listed in the previous section. Neglecting 
the effect of the inter-phase/interface, once we have decided the components, matrix and 
inclusion phases, inclusion concentration and their size, we can obtain different electrical 
properties for the mixture just by considering different geometries-distribution of 
inclusions ji^]. For this purpose, eleven different space-filling two-dimensional geometries 
are considered as presented in Fig. ^ The filler disks in two-dimensions are actually 
infinite unidirectional fiber inclusions in three-dimensions with the cross-sections in 
Fig. The structural information on the geometries are presented in Table [0 It is 
worth mentioning that some of the structures have similar limiting volume fractions q m , 



however, the coordination number Z 58( (number of closest neighbors) are different. The 



limiting concentrations are actually the percolation threshold I for a given 

regular geometry. Close to the percolation threshold local field is deformed from that 
of the average one (predicted by the effective medium approximation). The individual 
inclusions polarization is transformed from dipole approximation to multipoles. At the 
some concentration level the effective medium theories become unapplicable. 

If we return back to the two-dimensional structures by taking into consideration 
Fig.^and Tabled The geometries in Fig.^s; and^i have similar limiting concentrations 
(q m ~ 0.55) and coordination numbers. In such a case one can for example investigate 
the influence of the second nearest neighbors on the dielectric permittivity. Figs, 
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Figure 1. Two-dimensional structures considered in the simulations. The 
first three structures (a), (b) and (c) are the well-known square, triangular and 
hexagonal(honeycomb) networks, respectively. The others are generated from various 
space-filling tessellations. The regions marked with rectangles are the computation 
domains used in the finite element analysis. 



and ^ have similar limiting concentration (q m ~ 0.78) but coordination numbers. 
The geometries in FigGJ:, ^ and have on the other hand besides having similar 
limiting volume fractions (q m = 0.61), their coordination numbers differ from each 
other, introducing an interesting case to study the influence of structural differences and 
number of the closest neighbors on electrical properties. These significant differences are 
emphasized in the text by analyzing the electrical properties and the resulting spectral 
density functions of these structures with varied intrinsic properties of the constituents. 



3. Determining topology contributions 

3.1. Spectral (geometric) density function 

Bergman 0, 42, |4^{ has proposed a mathematical way for representing the effective 



dielectric permittivity e e of a binary mixture as a function of permittivities of its 
constituents, e\ and e 2 , and an integral equation, which includes the geometrical 
contributions. It is called the spectral density representation. The representation utilizes 
the Maxwell- Wagner- Sillars polarization and the geometrical contributions, which are 
calculated as a function of depolarization factors 

After the introduction of non-destructive measurement techniques and systems, 
such as electrical 110 or acoustic impedance S pect r o S co Py 0, the impedance of 



Table 1. Structural information regarding the considered two-dimensional geometries 
in Fig. n the coordination number Z, limiting volume fraction q m , the size of the 
computational geometry a x b in the maximum filler particle radius r m and the fraction 
of the filler in the computation domain n e . 



Fig. 1 


Z 






axb 




"e 


a 


4 


0.78 




lxl 




l 

4 


1) 


6 


0.90 




V3xl 




1 
2 


c 


3 


0.61 




\/3x3 




1 


d 


■5 


0.77 


1 ■ 


F V3xl 




3 
I 


e 


4 


0.68 




\/3x2 




3 
4 


f 


3 


0.45 


2- 


F V3x2 


F V3 


2 


g 


3 


0.55 


3- 


F\/3x2 


F2V3 


4i 


h 


3 


0.54 


2 


F a/2x2 


F V2 


2 


i 


5 


0.61 


1 ■ 


F V3x3 


F V3 


3 


j 


5 


0.84 


1 


F V3xl 


F V3 


2 


k 


4 


0.61 


1 + 


2V3x2 


F V3 


34 



materials (either pure or composite) could be recorded for various frequencies v. Then, 
the frequency could be used as a probe to obtain microstructural information with 
the application of the spectral density representation 



This can only be achieved if (i) no influence of v on the geometrical arrangement 
of phases is present [the geometry should be static at each frequency, meaning that 
no piezoelectricity exists in the constituents, and the elastic properties of the phases 
should be the same or similar to each other otherwise there would be a nonzero 
displacement vectors (deformation) in the composite], (u) the intrinsic properties of 



phases are known as a function of v. Numerical 69|, [7l|, l72[ and analytical |5 



53, 6 



approaches have been used and proposed to resolve the spectral density function for 



composites. Even a relation between t 
of relaxation times have been ascribed 



re spectral density function and distribution 



73l I74j . The distribution of relaxation times 



approach 0, 0| has previously been applied to illustrate the differences in random 
and regular structures 26 , • Although numerical approaches could be prefered over 
the analytical ones, which are empirical expressions and are not universal, they solve 
a nontrivial — ill-posed — inverse problem [72j. Here, the same numerical method as the 
distribution of relaxation times is applied to extract the the spectral density function of 
a binary mixture. The method is based on the Monte Carlo integration and constrained- 
least-squares algorithms. 

For a binary composite system with constituent permittivities S\ and E2, and 
concentrations q\ and g 2 , (qi + q 2 = 1), and with an effective permittivity e e , the 
spectral density representation is expressed as, 

A ei / Aji - Aj = I g J (x)[l + e- 1 A ji x]- 1 dx (2) 

where, = — £j, and is complex and frequency dependent. Aj is a constant, and 
depends on the concentration and structure of the composite. Several different notations 
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Figure 2. Parametric plot of the scaled mixture permittivity. The symbols are the 

analytical model of Maxwell-Garnett equation, and the solid lines ( ) are the values 

calculated from the spectral functions obtained from the proposed numerical method. 
The semi-circles from large to small corresponds to q 2 = {0.95, 0.7, 0.5, 0.3, 0.05}, 
respectively. The inset is the enlargement of the values close to the origin for 
q 2 = {0.05, 0.30}. 



69 1 . Here, we rearrange the 



have been used in the literature see Refs. |52[,|53[ and 
expression by Ref . , and obtain a similar one those used by Refs. 0] and [fifll ] . The 
spectral density function is g(x), and it is sought by the presented procedure. The 



spectral density function satisfies J gj(x)dx = qj 43, [7l| and J xgj(x)dx = qjqi/d 



where d is the dimension of the system. The shape of the inclusions in a matrix can 
also be related to d 
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631 ] . Finally, x is called the depolarization factor. 



The numerical procedure is briefly as follows: first the integral in Eq. (J2J) is 
written in a summation form over some number of randomly selected (known) re- 
values, x n G [0, 1]. This converts the non-linear problem in hand to a linear one with 
only gj n values being unknowns. Later, a constrained-least-squares is applied to get the 
corresponding g Jn -values: 

min 1 1 A - Kg in 1 1 2 and g jn > (3) 

where A is the left-hand-side of Eq. (j2J, and K is the kernel-matrix, [1 + e~ 1 AjiX n ]~ 1 . 
When this minimization is run over-and-over with new sets of values, most probable 
gj^-values are obtained. For a large number of minimization loop, actually the x-axis 
becomes continues — the Monte Carlo integration hypothesis. Finally, the weighted 
distribution of gj n versus x n leads g(x). This is achieved by dividing the x-axis in 
channels and averaging g Jn in each channel. 

Application of the numerical procedure to the Maxwell-Garnett expression [7^ 
should yield delta function distributions for g(x) The dielectric function for a d- 

dimensional (or composite with arbitrary shaped inclusions) Maxwell-Garnett composite 
is defined as 

s e = £l [l + dq 2 A 2 i (qi A 21 + de 1 )- 1 ]. (4) 
The resulting the spectral density function is then, 

gj (x) = 6[x-(l- qj )/d). (5) 
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Figure 3. Calculated spectral density distributions, which correspond to delta 
sequences. The spectral functions from left to right corresponds to q 2 = 
{0.95, 0.7, 0.5, 0.3, 0.05}, respectively. The corresponding (calculated) A 2 values 
are {0.002, 0.012, 0.029, 0.064, 0.358}, respectively, for the considered concentrations. 
The dashed lines (- -) show the positions of the actual delta-functions for the 
Maxwell-Garnett expression. 



We choose the following values for dielectric functions of the phases: E\ = l—i (IQQeqU))^ 1 
and e 2 = 10 — i {e^uS)" 1 . The left-hand-side of Eq. (J2J) without the constant A 2 is plotted 
for a 3-dimensional composite (d — 3 corresponding spherical inclusions) in Fig. |21 as a 
parametric plot of the imaginary part of A el /A 2 i against its real part. The graph is a 
semi-circle for the Maxwell-Garnett expression. In the figure, five different concentration 
levels are plotted, g 2 = {0.05, 0.3, 0.5, 0.7, 0.95}, the inset shows the enlargement close 
to the origin, which illustrates the low concentrations, g 2 = {0.05, 0.3}. The size of the 
semi-circles are proportional to the concentration of Phase 2. The analyses performed 
on the scaled effective permittivity, Eq. (j2J), with the help of the applied method yield 

the solid lines ( ) in the figure. 

The corresponding g(a;) are plotted in Fig. El on a log-log scale. In the figure, the 
expected locations of g(x) from Eq. (jSj) are also shown with dashed lines (- -). The 
g(a;)-distributions obtained are analyzed by the Levy distribution [76|, which generates a 
delta-sequence JE3]- The solid lines (— — ) illustrate the appropriate Levy distributions. 
Various parameters from the statistical analyses of the spectral density function, and 
their expected values are presented in Table El The concentration values, g^, calculated 
from the integration of g(x) without a-priori assumption are < 1% for the considered 
high concentrations, and it is around 5% for the lowest concentration, g 2 = 0.05. 
The localization parameter for the depolarization factor x, which is the most probable 
depolarization value, can be calculated by the integration of [1 — g(x)]/d or with the 
help of statistical analysis. The estimated depolarization factors x are within < 1% 
of the actual values stated by the proposed analytical expression ji^, 0] . Finally, the 
product of the concentrations gig 2 , the integration of 3xg(x), calculated have also very 
good agreement with those values expected from the definitions of the spectral density 
representation. 



Table 2. Comparison between the results of the proposed numerical approach and 
those of the Levy statistics and the given analytical the spectral density function for 
the Maxwell-Garnett effective permittivity expressions for various concentrations. The 
bars on the quantities indicate that they are calculated from the numerical results. 
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a Calculated using the resulting g2(x). Known from the definition of gj{x) — integral 

gj(x)dx is equal to this value. 
b The localization parameter for the calculated Levy distribution. The shape 
parameters and the amplitude of the Levy distributions are disregarded. 
c Known from the definition of the spectral density function for the Maxwell-Garnett 
expression, Eq. J5j|. 

d ^2-value calculated before the numerical procedure using Eq. 

Mean ^-value calculated during each Monte Carlo integration step in the numerical 
procedure, Eq. 

f Calculated using the resulting g2(x) and x- values. Known from the definition of 
gj(x) — the values is equal to the integral J Q 3xgj(x)dx. 
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Figure 4. Parametric (Cole-Cole) plots of dielectric permittivities with varied material 
parameters at different concentrations; match-composite at (a) q — 0.1 and (b) q = 0.4, 
and reciprocal-composite (c) q = 0.1 and (d) q = 0.4. The symbols represent the 
structures in Fig. as follows; (a:>, (b:o), (c:+), (d:<), (e:*), (£□), (g:v)> ( n:- ); 
(j:A) and (k:x). 

3.2. Numerical calculations 



In the simulations the geometries in Fig. ^ are used, and the electrical properties of the 
structures are calculated with the finite element method BB& The computational 
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Figure 5. Frequency independent parameters for the match composites for considered 
concentration levels (a) and (b) q = 0.1 and (c) and (d) q — 0.4 for the structures in 
Fig- ID ( a ) and (c) ohmic conductivity a versus permittivity at high frequencies e. (b) 
and (c) low frequency permittivity e + x(0) versus permittivity at high frequencies e. 



domains used in the simulations are marked in Fig. ^ with rectangles. These regions 
are the smallest repeating units (unit-cells) of the structures. The phase parameters are 
chosen to correspond to match- and reciprocal composites as presented in Table El For 
descriptions of match and reciprocal composites see Refs. 0, |3| . These parameters for 
the intrinsic electrical properties of phases yield relaxation times r around 1 s for the 
interfacial polarization in a dilute mixture. 

In Fig. |U dielectric permittivities calculated are plotted after the subtraction of the 
ohmic losses. At low concentrations q — 0.1, it is not important if the composite is match 
(Fig. Ht) or reciprocal (Fig. Ht), dielectric dispersions of structures are almost the same 
except the structure with the lowest limiting concentration, Fig. Qf. Because of the low 
concentration, the polarization of each fiber is not affected by the neighboring fibers, so 
the resulting polarization of the whole structure/composite material is as if each fiber is 
alone, and does not influenced by the nearest neighbors polarization. However, for the 
structure in Fig. the influence of the neighboring fibers are pronounceable. At high 
concentrations of the fibers (phase 2), the differences between each structure start to be 



Table 3. Material parameters adopted in the finite element simulations. 
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Figure 6. Frequency independent parameters for the reciprocal composites for 
considered concentration levels (a) and (b) q = 0.1 and (c) and (d) q = 0.4 for the 
structures in Fig. ^ ( a ) an( i ( c ) ohmic conductivity a versus permittivity at high 
frequencies e. (b) and (c) low frequency permittivity e + %(0) versus permittivity at 
high frequencies e. 



significant due to influence of the polarization of the neighbors in Fig.EJ} andEJi. It is not 
obvious to the naked-eye that all the calculated dispersions are deformed Debye 78| . It 
has previously been shown that for two-dimensional composites jjfiiEsI and in dispersive 
systems (z^l non-Debye relaxations would occur. Therefore, the dispersion data are 
analyzed by the Havriliak-Negami 80] empirical formula. 

e hn (oj) = e + X (0)[l + (uuTx) a ]- p + o{ie uu)- 1 (6) 

Here, a and (3 are fitting parameters that define the shape of the dispersion, r and x(0) 
are the position (relaxation time) and its amplitude (dielectric strength), respectively. 
This equation is the same as Eq. x{ u ) being the Havriliak-Negami empirical formula. 
In Tables 0] and Havriliak-Negami fits are presented together with the fitting errors. 
The error is calculated from the following expression, 

error = ^{[3^ e - e hn )/K(e e ) 2 } + [3(e e - e hn )/Z(e e ) 2 }} (7) 

The calculated errors indicate the fitness of the model function in Eq. (jUJ). Only at high 
concentrations, the structures with low limiting concentrations q m or in other words 
low percolation threshold show poor fitting results, which indicate that actually the 
Havriliak-Negami equation is not suitable for such an analysis. Although some of the 
(3 values are over 1, the dispersion satisfies the condition a < 1 and a(3 < 1 j^- 
The relaxation times r obtained from the curve fittings for match composites at low 
concentrations are approximately the same. The shape of the relaxations are also 
Debye-like. However as the concentration is increased due to the structural differences 



Table 4. Havriliak-Negami curve fitting parameters for the considered match 
composite structures. The last column is the error calculated using Eq0 
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and the interaction between inclusions, each geometry has a significant relaxation time 
t. Moreover, the shape parameters of the Havriliak-Negami empirical formula changes. 
As we consider the fitting results of the reciprocal composites, they illustrate significant 
differences between the structures. Keep in mind that the concentration of the inclusions 
are the same in Tables |U and the span of the relaxation times r and x(0) show 
the influence of the structure. If we have the data as a black box, it would not be 
possible to calculate the concentration of the each phase as presented in the previous 
section, however, if we instead use the spectral density representation we would obtain 
valuable information regarding the structure of composite as presented in Table El The 
relaxation times of the match and reciprocal composites are altered and they are faster 
for the reciprocal composites due to more conductive matrix phase. Finally, the main 
difference between the reciprocal and match composites are the size of the dielectric 
strength which is higher for the former case because of the Maxwell- Wagner-Sillars 
polarization. 

In Fig.Elandinifrequency independent parameters of the composites are presented to 



Table 5. Havriliak-Negami curve fitting parameters for the considered reciprocal 
composite structures. The last column is the chi-square error calculated using Eq0 
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give some more information on the structures and to be used to identify the dispersions 
in Fig. EJ For match composites at low concentrations there are small dissimilarities, 
illustrated in Fig. and Et>- As presented previously when the concentration of the 
inclusions in the composite systems get close to the percolation threshold the structural 
differences dominate the electrical properties. This is shown in Fig. Efc and Eli, in which 
structures in Fig.^F,^; and 0i separate from the others. When the conductivities of the 
phases are interchanged and a reciprocal composite is formed even at low concentrations. 
The discrepancy between different structures are significant this is due to the increase 
of polarized charge at the interface. 

Once the dielectric dispersions are known together with the intrinsic dielectric 
properties of the constituents, the data can be converted to the spectral density 
representation form, the scaled dielectric permittivities A^- in Eq. (j2J . This is illustrated 
in Fig. [7| It is graphical to see the structural differences in the scaled permittivities. 
The smaller semi-circles in the positive quadrant are for match composites for q = 0.4. 
The negative imaginary parts are obtained for reciprocal composites for the same 
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Figure 7. Parametric representation of the scaled composite permittivity for selected 
three structure from Fig. ^ The negative imaginary parts are obtained for reciprocal 
composites. 

concentration. It is significant that Fig.^i and^; change behavior as the real part of the 
scaled permittivity increases. Moreover, it is similar to Fig. that the structure with 
the lowest q m has the largest deviation of all, one very significant feature of this structure 
is that the flat part marked with the solid line (— — ) in the figure, which illustrates the 
deviation from the simple Debye relaxation. As a note, the Debye relaxation observed in 
layered-structures correspond to a pole at x = 1 or x = in the spectral representation, 
which is actually a delta function S(x — 1) or 5(x) depending on the direction of the 
applied field and the layered structure. 

The analysis of the scaled permittivities with the help of the developed numerical 
method resulted in the spectral density functions presented in Fig. |H| It is clear that the 
structure in Fig. ^ is different than the others at low concentrations, Fig. |H^l, because 
its peak is not located near the solid line (the expected position for the depolarization 
factor from the Maxwell- Garnett expression). Although it is previously stated many 
times in the present text, all the other structures yield the same spectral distributions, 
emphasizing that they are not actually different. This is also a verification of the 
effective medium approach. (A similar result can also be drawn from Fig. andEfc.) 
When the high concentration case (q = 0.4) is taken into consideration, Fig. |Ht>, there 
are significant differences between the structures. Now only geometries Fig. ^p, 
and C] yield similar spectral density distributions, which coincide with the suggested 
depolarization factor of the effective medium theory [solid line (— — ) and Eq. 
This reveals that these structures are in fact not close to the percolation threshold. The 
other structures with lower percolation thresholds (q m < 0.8) spread in the spectrum. 
The structure in Fig. approaches to zero depolarization factor, giving a sign for start 
of a percolation-touching of fibers to form a continues path from one side to other. In 
Table® calculated positions of the structural parameters are presented. There are little 
discrepancies between values obtained whether the composite is match or reciprocal. 
The table also verifies the statements regarding position of the depolarization factors, 
x and the percolation thresholds. 



Table 6. Peak positions of depolarization factors 3? calculated with the Levy statistics 
for the geometries in Fig. ^ The numerical method to extract the spectral density 
function is applied both to match and reciprocal composite configurations. The 
expected values for depolarization factor from the effective medium theory for q = 0.1 
and q = 0.4 are 0.45 and 0.3 [calculated from Eq. (0 with d = 2], respectively. 
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4. Discussions 

In this section, I would like to take three examples and with the help of the information 
presented in the previous section characterize the differences between structures. For 
this analysis as described in § 12.21 the structures in Figs. and [T]i, Figs[T|s; and HJi, 
and FigsHfc, d and are selected. 

First, let us consider the structures in Figs. EJi and Qi, which have nearly the 
same percolation threshold but the number of the nearest neighbors, Table H Since 
the considered concentration values are away from the percolation thresholds, the 
Havriliak-Negami empirical formula and the spectral representation approach resulted 
in a very similar way for the both concentration cases for the match composites and 
low concentrations of reciprocal composite arrangements. For high concentrations 
q = 0.4, the reciprocal composites of these structures have slightly different curve fitting 
parameters as presented in Table and the errors in the fitting are higher if compared 
to the ones at low concentrations, Tables 0] and EJ The depolarization factors from 
the spectral density analyses in Fig. |H1 and Table |B1 on the other hand illustrates minor 
differences for both match and reciprocal cases. The spectral representation states that 
the obtained spectrum contains contribution of the pure geometry, it is clear that we 
are not able to see the difference between these two structure when we compare the 
obtained depolarization functions, Table |H1 

Secondly, we concentrate on the structures in Figs, ^s; and QJi, which have similar 
percolation threshold and the same number of the nearest neighbors. For both 
simulation considerations, match and reciprocal, at low concentrations these structures 
yield similar dielectric dispersions, see Tables 0] and El However, at high concentration 




Figure 8. Calculated spectral density distributions, (a) At low concentrations 
density distributions are similar and like delta sequences, (b) At high concentrations, 
distributions are broader. The solid lines are the expected positions for the spectral 
density distributions from the effective medium theory Eq. Ip|l. The solid and dotted 
lines are calculations from match and reciprocal composites, respectively. 



of inclusions q = 0.4, only the relaxation shape parameters a and (3 and the position 
r of the dispersions are alike in the curve curve fitting analysis with Eq. (jSJ. The 
dispersions amplitude x(0)j permittivity at high frequencies e and ohmic conductivity 
a differ from each other. For reciprocal composite case the calculated error is large, 
however, the same is not true for the match composite, therefore the Havriliak-Negami 
empirical formula can be used to differentiate differences between structures. However, 
the spectral representation approach as presented in TableElyield slightly different values 
for the depolarization factors, which indicate that the structures are separable. We can 
also speculate that the differences in the frequency independent parameters might have 
been due to the influence of the second nearest neighbors. One other explanation can 
be the long range interaction of the clusters of inclusions, this can be seen as the change 
in the orientation of inclusions in groups of four in Figs. ^ and QJi- 

Finally, the three structures in Figs. ^ and ^ and ^ have the same percolation 
threshold (q = 0.61) but different coordination numbers, 3, 5 and 4 respectively 
(Table HJ. At low concentrations far away from q m , these composites yield the same 



Table 7. A summary of the dielectric relaxation in binary composite depending on 
the intrinsic electrical properties of constituents and topology. The inclusions in the 
simulations on regular structures have been mono-dispersed as in the present study. 
The deduction are from the simulations performed by the author. The concentration 
low and high means relation to the limiting concentration q m in regular topologies and 
percolation threshold in the disordered ones. S-reciprocal is the slightly reciprocal case, 
and E-reciprocal is the extreme reciprocal case. For E-reciprocal composites, dielectric 
dispersion of the system exhibits in low-frequency-dispersion (lfd) or in other words 
constant phase angle behavior [2^ . 
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spectral density function. The application of the Havriliak-Negami function also 
confirms the findings. At high concentration, q = 0.4, the differences between the 
structure become visible in both analyses, the spectral density representation and the 
Havriliak-Negami curve-fitting. The match and reciprocal composite cases indicate that 
the structure with the largest coordination number (Fig.^) yield low permittivity values 
compared to the others, Tables |U and El When the conductivity is considered, the same 
is observed for the match composite, Table 01 which is expected due to the similarity of 
the capacitive and resistive problem. The reciprocal composite of the same structure on 
the contrary results the opposite and yield a highest conductivity o of all three, Table El 
The two structures Figs. ^ and ^ have similar electrical properties, when the error 
in the curve-fitting are taken into consideration in Tables |U and the differences are 
within the error. The curve fittings reveal that the dispersions obtained are non-Debye 
type and depending on the intrinsic properties of the phases can be symmetrical or 
non-symmetrical. 

As a remark, the dielectric dispersion in two-dimensional binary composites can 
be summarized as in Table The items in the table have been obtained purely from 
numerical simulations, published previously elsewhere 0, 13 26, 28, 30, At low 
concentrations of inclusions the dielectric dispersion are in the Debye form, and they 
can be analyzed by the effective medium theories for a given shape of inclusion, when 
the inclusions are mono-dispersed. The term "low concentration" in Table [7| is meant 
that the concentration of inclusions is relative to the limiting concentration, the level 
at which the inclusions start to touch each other and form a continues chain from 
one-side of the composite to the other one. In the disordered/random structures, the 
low concentration cases would then be assigned relative to the percolation threshold. 



In disordered/random composites the dielectric dispersions are non-Debye type, and 
depending on the intrinsic electrical properties of the constituents even low frequency 
dispersion (lfd) or also known as constant phase angle [l|, 0, EH 0, EE Edl EB type 
of dispersions are obtained on a wide-frequency range. The frequency range of the lfd 
dispersions can be increased if the number of components (inclusion phase) or the largest 



to smallest length scale ratio are increased in the numerical simulations 24] . Although, 
only the position of the depolarization factors are presented and considered in the data 
analysis and discussions, the shape of the spectral functions also provide information 
regarding the micro-structure of the composite. As stated by Lysne 73j . the position 
of the depolarization factor in the spectral density representation is correlated with the 
relaxation time of the dielectric dispersion. Therefore, one can in principle use the 
distribution of relaxation times , however, it would be cumbersome, since there 

is no trivial quantity as the depolarization factor in the spectral density representation 
to correlate the relaxation times. 



5. Conclusions 



In this paper, an insight to the electrical properties of binary composites are given 
with the help of numerical simulations. It is presented that if we are able to control 
the manufacturing process in preparing materials for desired electrical properties, the 
computer simulation help us not only to characterize but to design materials as well. 
Eleven selected space filling geometries (lattices) are used in the numerical simulations. 
The structures are "ideal composites", in which we can change the intrinsic properties 
and distribution of inclusions. The dielectric response for each structure is calculated for 
two concentration levels of inclusions. The obtained dielectric dispersions are analyzed 
by the Havriliak-Negami empirical formula and the spectral density representation 
method. It is concluded that for low concentrations of inclusions (lower than the 
limiting concentration or the percolation threshold for a given system), the composite 
systems exhibit Debye-type dispersions, which can be modeled by empirical formulas for 
mixtures. At concentrations close to the percolation threshold, the structures illustrate 
deviations from the Debye-type relaxation. The spectral density approach illustrate 
some advantages over the effective medium approach and using empirical formulas 
for the dielectric dispersions obtained when the systems have high concentrations of 
inclusions. It is now time to create look-up tables for high concentration composites 
with known micro-structures in order to compare experimental and computer simulated 
dielectric dispersion of composites. With the help of the look-up tables, we would be able 
to know the internal structure, such as the distribution of inclusions. Thus, the dielectric 
spectroscopy data on composites together with the spectral density representation can 
be used as a probe to provide quantitative information on the micro-structure of the 
composite. 

As a scheme or a next step, it is valuable to prepare composites with known internal 
structures and verify experimental results with numerical simulations or the opposite. 



With the improving material preparation techniques, this should be a task for both 
theorists and experimentalists. 
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